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INTRODUCTION 


This semi-annual status report covers the period from June 1, 1977 
to November 30, 1977 and contains a review of the research and applica- 
tions, completed or in progress, as funded by the Office of University 
Affairs, NASA and conducted by Purdue University, Laboratory for Appli- 
cations of Remote Sensing. 

This reporting period marks the first half of the fifth year of 
funding for a proposal entitled "The Applications of Remote Sensing 
Technology to the Solution of Problems in the Management of Resources 
in Indiana." As indicated in this title, the purpose of this work is 
to introduce remote sensing into the user community within the state of 
Indiana. The user coimiunity includes those local, regional and state 
agencies involved in the decision monitoring and/or managing processes 
of the state's resources. 

In order to carry out this work it is not only necessary to initiate 
projects with these agencies but also it is necessary to meet with and 
provide information to as many people and groups as well as agencies as 
possible. During the past six months numerous meetings are held with 
many different groups. 

Among the groups that were contracted and received information 
about th-is program were: 

Area Planning Commission, Tippecanoe County 
Area Planning Commission, Boone County 
St. Joseph County Area Plan Commission 
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Michiana Council of Governments 
Indiana Geological Survey 
U.S. Forest Service 

Tipton County Commissioners and Engineers 
Indiana Department of Natural Resources 

a) Division of Reclamation 

b) Division of Forestry 

c) Division of Properties, Fish and Wildlife 

d) Soil and Water Conservation Conmittee 

Soil Conservation Service. 

Listed below are the projects that are reported in this document; 
Soils Inventory 

Forestry Demonstration Project 
Coastal Zone Management 
Tippecanoe County, Land Use Project 
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SOIL INVENTORY PROJECT 


INTRODUCTION 

The acceleration of the National Soil Survey program and the 
production of useful, high quality soil surveys in Indiana are among 
the prime goals of the USDA/Soil Conservation Service and the Indiana 
Department of Natural Resources Soil and Water Conservation Committee. 

The wide use of soil surveys for engineering and planning purposes in 
addition to agricultural uses has resulted in many specific questions 
concerning the physical nature of the different soil units depicted on 
soil maps. In order to provide the details necessary to understand the 
landscape composition and to provide interpretation of soil maps for 
specific uses, information of a quantitative nature is needed. To ac- 
complish this task all avenues are being considered, including remote 
sensing technology which can provide quantitative measurements through 
computer analysis of Landsat multispectral scanner (MSS) data. 

OBJECTIVE 

The overall objective of this task is to determine the applicability 
of using computer analysis of Landsat multispectral scanner data in 
accelerating and improving the quality of the soil survey program in 
Indiana. 

To evaluate the usefulness of the data the following specific 
studies were initiated: 

1. Evaluate the usefulness of spectral soil maps produced from 
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multispectral scanner data using pattern recognition 
techniques as quality control in soil surveys and as 
a means to evaluate quantitatively the soil mapping 
unit composition. 

2. Investigate the possibility of producing high quality 
general soil maps using false color Landsat imagery as 
the base map. 

3. Develop a soil parent material map using multispectral 
resource data. 

4. Determine the feasibility of producing a spectral soil 

map on a county-wide basis with its accompanying manuscript 
and evaluate the utility of this type of soil survey report 
to user groups. 

5. Evaluate the usefulness of superimposing computer classi- 
fication results upon aerial photobase maps in order to 
gain the benefit of the landscape prospective, 

STATUS 

Research in this area is now being conducted solely in Jasper County, 
Indiana. As stated in past reports a data set consisting of Landsat MSS 
data geometrically corrected and precision registered to ground center 
points at a scale of 1:15840 has been prepared. This data has been over- 
laid with a set of boundaries delineating the four types of parent materials 
found within the county: 

1. Glacial till 
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2 . Outwash 

3. Outwash over glacial till 

4. Lacustrine 

These boundaries were developed by visual interpretations of 
simulated false color imagery and single band grey scale images re- 
presentations of the Landsat data. 

Twelve quarter sections were randomly selected in three different 
parent materials, namely outwash, lacustrine and glacial till. These 
quarter section plots have been field mapped of a scale of 1:7920. At 
this scale, contrasting detail, down to one hectare in size has been 
delineated. These sites have been located on the Landsat data, and they 
will be used for correlating spectral data with soils of known parent 
material and drainage characteristics. If necessary, addition quarter 
section areas will be mapped to assure representation of all the soils 
in the county. 

Four spectral classification of Jasper county have been completed. 
All four techniques used an unsupervised approach however these classi- 
fications differ by the manner in which the training statistics were 
developed. The four methods of selecting training statistics are 
listed below: 

1. a) Subjective selection of areas representing the 
spectral variability within each parent material 
area; 

b) aggregation of statistics into one set of statistics 
for the entire county; 



6 


c) classification of entire county. 

2. a) Subjective selection of areas representing the 

spectral variability within each parent material 
area; 

b) development of individual training statistics 
for each parent material area; 

c) classification of each parent material area 
using its representative statistics. 

3. a) Objective sampling of the entire county using 

.every eleventh line and eleventh column of data; 

b) development of training statistics for the entire 
county; 

c) classification of the entire county. 

4. a) Objective sampling of each parent material area 

using every fifth line and fifth column of data; 

b) development of training statistics for each parent 
material area; 

c) classification of each parent material area using 
its representative statistics. 

Evaluation of the four classifications are currently being conducted 
by comparing them to the twelve quarter sections previously mentioned. 

The two classifications best representing the soils within the county will 
be taken to the field for further evaluation. Selection of the final 
classification will be made after the field evaluation and correlation of 
spectral classes to soil classes has been completed. 
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The final product for Jasper County will consist of a digitized 
multi spectral soil map produced at a scale of 1; 15,840 printed on clear 
mylar. This will underlay a matte halfone film positive of the aerial 
photo field sheet. In this manner the multispectral map will be used 
to guide the soil scientist In determining where to place soil boundaries 
for the various mapping units and to also point out specific locations, 
in the landscape, that are different and where additional investigative 
borings need to be made. 

A more detailed report on the results of this effort will be possible 
at the end of the next reporting period at which time the field correlation 
phase will have been completed. 
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TIPPECANOE COUNTY LANDUSE STUDY 


INTRODUCTION 

Great need exists to inventory landuse on the county level for 
planning purposes. As urban areas continue to grow and engulf rural 
settings it becomes imperative to inventory land cover and use all 
available information in landuse planning of rural lands. 

The Tippecanoe County Area Planning Commission (APC) has expressed 
a growing interest in utilizing remote sensing data as a planning tool 
along with its standard maps and black and white aerial photography. 

As a consequence a demonstration project was devised to see what aid 
digitally derived landuse maps might benefit the APC. 

STUDY SITE 

Tippecanoe County is located in the west-central part of Indiana 
(see figure 1). It is a rectangle 501 square miles, or 320,640 acres, 
in area. The county is approximately 7055 agriculture, 6 % urban and 24% 
other land uses. 

Numerous housing developments are scattered over the outlying 
rural areas and strip development along major highways abound. 

OBJECTIVE 

The major objective of this demonstration project was to utilize 
multispectral scanner data (MSS) with its .5 hectare (1.1 acre) reso- 
lution to categorize each .5 ha cell as to its land use with particular 
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emphasis in the rural areas. Project personnel who were experienced in 
previous landuse studies also were aware that rural-urban interface 
areas were unique problem areas and would need special attention during 
the study. 



Figure 1. Location of Tippecanoe County, Indiana which was used in 
this landuse danonstration study. 
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APPROACH 

The overall approach to this study was to produce spectral 

analysis of Tippecanoe County, Indiana utilizing Landsat MSS data and 

% 

computer implemented pattern recognition analysis techniques. These 
results were evaluated for accuracy by means of representative test 
fields. 

In addition, a secondary approach was to delineate urban areas in 
the county and digitize their boundaries were to be used as an additional 
channel of data. These urban boundaries were to be utilized along with 
the spectral data to allow for classification of the urban protions of 
the county using the urban class statistics and the rural areas using the 
rural class statistics. These analysis results were to be evaluated by 
the same test fields used in the earlier evaluation. The two classifi- 
cation results were to be compared and differences noted and interpreted 
where possible. 

ANALYSIS PROCEDURE 

A data set over Tippecanoe County for June 1976 (.Landsat scene i.d. 
2515-15411) was selected because of the in-house availability of the 
data and the fact that March 1976 1:24,000 scale black and white aerial 
photographs were available through the Tippecanoe County Area Planning 
Commission. These MSS data were geometrically corrected and registered 
to portions of 16 U.S.G.S. 7% minute topographic maps of the county at 
a 1:24,000 scale. 

Representative areas containing the different covertypes found in 
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the county were handpicked from a digital image displayed on a black 
and white CRT screen with the aid of the 1:24,000 aerial photographs. 
Sixteen initial training areas were processed using the LARSYS CLUSTER 
algorithm to determine the natural spectral structure within each area. 
Table 1 lists these areas and their line and column coordinate location. 
Figure 2 illustrates a map showing the location of the training areas. 


Table 1. List of the initial 16 training areas with their coordinate 
location and major landuse group assignment. 


Training 

Area 

Area 

Training 

First 

Last 

First 

Last 

Number 

Type 

Area 

Line 

Line 

Column 

Column 

1 

Rural 

HERSHEY 

213 

248 

497 

545 

2 

Urban 

ORCHARD 

288 

314 

409 

437 

3 

Urban 

MALL 

327 

349 

377 

401 

4 

Urban 

DWNTWN 

285 

308 

297 

330 

5 

Urban 

OLD-RES 

304 

324 

336 

357 

6 

Urban 

NEW-RES 

339 

355 

333 

361 

7 

Rural 

MCCUTCH 

421 

468 

332 

387 

8 

Rural 

WST PT 

403 

433 

113 

150 

9 

Rural 

Trailer 

351 

380 

221 

258 

10 

Rural/Urban 

LILEY 

329 ■ 

344 

238 

285 

11 

Rural 

DAYTON E 

356 

387 

540 

566 

12 

Rural 

DAYTON 

348 

379 

484 

513 

13 

Rural /Urban 

WATER/WD 

280 

310 

126 

156 

14 

Urban 

SUBD 

256 

266 

156 

171 

15 

Urban 

IND/HWY 

223 

231 

150 

160 

16 

Rural 

FLDPLAIN 

307 

329 

189 

207 


The spectral cluster classes for each training area were related 
to specific group covertypes with the aid of the black and white aerial 
photographs and the U.S.G.S. topographic maps. Training areas were then 
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Figure 2 _ Map of Tippecanoe County showing location of the 16 training areas 
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assigned to a major landuse group (urban, rural, urban/rural) and their 
respective cluster classes were grouped together into single sets of 
urban classes, rural classes, and urban/rural classes. 

Each of the resultant set of classes was tested for separability 
between classes within the set by utilizing the SEPARABILITY processor. 
The outcome of this step was a recommended set of 20 urban classes, 12 
rural classes and 7 urban/rural classes as listed in Table 2, 


Table 2. List of the 39 spectral cluster classes for the urban, 
rural and urban/rural class groups. 


Rural Classes 


Agriculture 1 

Bare Soil 1 

WetSoil 1 

Agriculture 4 

Agriculture 2 

Light Soil 

Tree 1 

Pasture 1 

Agriculture 3 

WetSoil 2 

Saturated Soil 1 

Bare Soil 2 

Urban Classes 

Commercial 1 

New Residential 

Sparse Vegetation 

Residential 2 

Trees 2 

Trees 3 

Wet Vegetation 

Mall 

Wet Soil 3 

Wet Agriculture 1 

Saturated Soil 2 

Commercial 4 

Comnercial 2 

Shallow Water 

Commercial 3 

Water 2 

Residential 1 

Wet Soil 4 

Street 

Water 1 

Rural /Urban 

Light Soil 2 

Wet Agriculture 2 

Sparse Agriculture 2 

Water 3 

Pasture 2 

Trees 4 

Bare Soil 3 

Industry 


These 39 classes were than grouped together to a single statistics 
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set and the separability between these new classes were calculated. 

Those classes which were similar spectrally were merged together and 
a final set of classes for the whole county was obtained. This final 
list of classes is given in Table 3. 

The statistics (means and covariances for each class) for these 
final classes were then used along with the CLASSIFY-POINTS processor 
to classify every .5 hectare cell (data point) in Tippecanoe County. 

The results map is shown in figure 3 and acreage tables given in Table 4. 


Table 3. List of final 20 classes used to classify 
Tippecanoe County. 


Gravel 
Bare Soil 
Light Soil 
Pasture 

Wet Agriculture 1 
Wet Agriculture 2 
Saturated Dark Soil 
Sparse Vegetation 
Agriculture 
Trees 1 


Trees 2 
Residential 
New Residential 
Streets 
Commercial 1 
Commercial 2 
Industry 
Shallow Water 
Water 1 
Water 2 


SPECTRAL CLASSIFICATION EVALUATION 

Preliminary random test fields were selected from all parts of 
Tippecanoe County. These fields were photo- interpreted to select those 
which contained homogenous identifiable covertypes. The selected fields 
used for testing classification results contained 1138 sample data points 
which yielded an overall classification performance figure of 86.6%. 

Table 5 lists the class groupings that were tested (agriculture, forest, 
urban) and their respective performances. 


ORIGINAL PAGE IS 
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Figure 3. Classification results maps of 
Tippecanoe County. Indiana. 




Well drained agricultural land 

Pasture land 

Tree covered areas 

Poorly drained agricultural land 

Residential areas 

Industry/Coflinercial areas 

Water 



Table 4. Spectral classification results of Tippecanoe 
County, Indiana, in tabular format. 


Class 

Total Acres 
Classified 

Percent 
of County 

Gravel 

4572.15 

1.4 

BSOIL 

36900.3 

11.1 

LTSOIL 

1484.04 

0.4 

Pasture 

27736.39 

8.3 

WETAG 1 

36834.53 

11.1 

WETAG 2 

52862.43 

15.8 

SATSOIL 

44844.44 

13.5 

SPRSVE6 

7932.60 

2.4 

Ag 

43738.91 

13.1 

Trees 1 

32223.14 

9.7 

Trees 2 

19620.31 

5.9 

Resident 

13591.81 

4.1 

New Res 

586.23 

0.2 

STREET 

4292.88 

1.3 , 

Comerc 1 

2436.09 

0.7 

Comerc 2 

431.60 

0.1 

Industry 

1332.87 

0.4 

SHWATER 

544.69 

0.2 

Water 1 

115.4 

0.0 

Water 2 

963.59 

0.3 


333044.4“ 100.0 


DELINEATION OF URBAN AREAS 

The second step in this study involved delineating the urban areas 
throughout the county. This was accomplished by examining ARC zoning 
maps along with black and white aerial photographs. All cities, in- 
corporated towns and major housing subdicirions were outlined on 7% 
minute topographic maps so that they might later be digitized and their 
coordinates recorded on magnetic tape. A total of 13 cities and towns 
and 15 subdivisions were delineated for this study. A classification 
of the county using this boundary information has not been completed, 
but will be presented in the next semi-annual report. 
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Table 5. Classification Performance. 


TEST CLASS PERFORMANCE 


NUMBER OF SAMPLES CLASSIFIED INTO 



NO OF 

PCT. 

GROUP 

SAMPS 

CORCT 

1 AG 

819 

93.7 

2 FOREST 

148 

86.5 

3 URBAN 

171 

52.6 

TOTAL 

1138 



AG 

FOREST 

URBAN 

WATER 

767 

31 

18 

3 

18 

128 

. 2 

0 

64 

14 

90 

3 

849 

173 

no 

6 


OVERALL PERFORMANCE ( 985/ 1138) = 86.6 

AVERAGE PERFORMANCE BY CLASS (232.8/ 3) = 77.6 



18 


CONCLUSIONS 

The completion of the first part of this landuse study has made 
possible the following conclusions: 

* Classification performance of the spectral analysis 
was significantly higher than was expected consider- 
ing the scattered housing and strip developments. 

* Numerous spectral classes which exist in both urban 
and rural environments in reality have different 
land uses (ex. park grassland - pasture, disturbed 
soil in housing developments - plowed farm soils, 
city parks - rural tree lots, etc.). 

* Urban-rural interface zones are complex areas need- 
ing concentrated attention during analysis procedures. 

* The addition of digitized urban boundaries should 
greatly aid discrimination between spectrally 
similar urban and rural classes. 

The protion of this study involving utilization of digitized urban 
boundaries appears to be an exciting and promising step toward develop- 
ing a usable landuse planning tool. 
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AN EVALUATION OF INDIANA'S COASTAL 
ZONE FOREST RESOURCES 
A Study of the Impact of Remote Sensing 
to Forest Inventory 


INTRODUCTION 

This report describes the results from a forest inventory con- 
ducted within the boundary of the Indiana Coastal Zone Management 
Area (fig. 1). The materials documented here are from the second 
phase of a two-year study which began in 1976. This inventory 
utilized combined inputs from remote sensing and ground based tech- 
nologies to characterize the extent and quality of the forest re- 
source. These results represent a quantative frame work upon which 
to build forest management policy in the Coastal Zone Management Area. 
Furthermore, this study demonstrates the feasibility of using Landsat 
collected data in conjunction with on-going survey to provide im- 
portant and necessary Information to direct the development of Indiana's 
forest resource. 

Background 

In response to Public Law 92-583 (the Coastal Zone Management Act 
of 1972) personnel from the Indiana Department of Natural Resources' 
(IDNR) Division of Forestry, and the Laboratory for Applications of Re- 
mote Sensing (LARS) cooperated in a forest inventory of the Coastal 
Zone Management Area (CZM). LARS provided maps and statistical data for 
forest cover through machine-assisted analysis of Landsat multispectral 
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Figure 1. Location of the Coastal Zone Management area in northern Indiana. 
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scanner data. Using this input IDNR inventory the timber lands for 
the purpose of categorizing the forest resource. Results from this 
first phase of the study were published in "Technical Report No. 101." 

The second phase of the CZM forest inventory was initiated in the 
spring of 1977. The intent of this new phase was to provide a more 
detailed quantification of the forest resource base. Specifically, 
information on the potential and present productivity of forest lands, 
and a measure of their susceptibility to erosion was to be provided. 
Given these inputs, broad management objective have been defined, along 
with unique areas which require special management consideration. 

Objectives 

The objective of Phase II was: 

To further amplify, and quantify forest 
resource information collected during 
the 1976 inventory. 

In order to accomplish this objective, we needed to: 

• Identify the potential productivity, of 
forests in terms of board foot growth/ 
acre/year, based on soil productivity. 

« Determine the potential erodability, of 
forest soils. 

• Classify the existing resource in terms 
of annual production, emphasizing areas 
of high commercial volume productivity, 
recreation potential or other unique 
attributes. 



22 


These objectives address specific requirements which were necessary 
to formulate management directions for the CZM. The scope of this 
inventory did not allow for the intensive collection of data neces- 
sary to develop managanent plans for a small ownership. Therefore, 
any inferences made from these data to individual forests tracts would 
be subject to excessive sampling error. Individual forest owners who 
desire management plans for their land should consult with the INDR - 
Division of Forestry. 

APPROACH 

The 1977 survey utilized the previously classified Landsat data, 
and current soil maps. The forest soils were identified from the Land- 
sat data and were digitized in 2.5 acre cells by section and township. 
Soil Conservation Service information on woodland suitability and land 
capability forms the bulk of a data base created at LARS. Field data 
to supplement the 1976 survey was also collected. IDNR processed the 
field inventory data and LARS manipulated the data base. Both organi- 
zations worked together to complete the inventory and develop recommen- 
dation. 

Methods and Materials 

Landsat data classified by LARS Staff early during 1976 (fig. 2) 
was the basis for the forest inventory. Fundamentally, the Landsat 
classifications were used to locate forest land. Four classes, based 
on distinct spectral response characteristics were identified. Forest 
inventory data were sampled based on a proportional allocation of acreage 
by spectral class. These survey data were reduced by IDNR, and the 
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results analyzed to provide a profile of the CZM forest resource. 

An important facet of this study involved combining the results 
of the Landsat classification with soils data. The soil data was 
derived from current soil maps for townships within the CZM. The 
soils information from the map can be related to indicies of potential 
woodland productivity, and potential erodibility. The Landsat and 
soils data were combined in a grid-cell data base, which conceptually 
illustrated in figure 3. 

The data base was created by encoding the following information 
about a 2.5 acre cell into a digital computer record: 

County Name 
Township 
Range 
Section 

Row Number (1-16) 

Column Number (1-16) 

Spectral Class of Forest 
Soil Symbol 
Slop Class 
Erosion Class. 

A grid, figure 4, at two scales, 1:15840 and 1:24,000 was used 
to compensate for the difference in scales between the Landsat and 
soils maps. 

Line printer classification maps for the forest classes only (fig 
5a) were manually annotated to show township, range, and section lines 
The 1:24,000 grid was used to locate the coordinates of forest tracts 



25 




Soils 



Forest classification 


Figure 3. Conceptual data base showing political, soils and 


vegetation layers. 





1 MILE 

@ 1:15,840, 1 mile = 4 inches - 

0 1:24,000, 1 mile = 2.64 inches 
ISECTION 




Figure 4. Example of grid use to digitize information 


about soil and vegetative cover for each 
section of land with the CZM. 
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consisting of a mininiuin of 2 pixels (about 2.3 acres). Soils infor- 
mation for these coordinates were than read from the soils map (fig. 
5b) and encoded to a computer file, (fig. 5c and Appendix V-1) and ex- 
ample below. 


FILE. . 

. THNSHP 

373 

PI 


LAPORTE 

37N03W030101 

3 


MX 

LAPQRTE 

37N03W030201 

2 


BTA 

LAPORTE 

37N03M03G505 

3 


MRB2 

LAPORTE 

37N03W03C605 

3 


MR02 

LAPORTE 

37M03W030701 

4 


CHB 

LAPORTE 

37N03W03O7O5 

3 


WT 

LAPORTE 

37N03W030901 

3 


BTA 

LAPORTE 

37N03W030903 

3 


RR 

LAPORTE 

37N03W031001 

4 


FH 

LAPORTE 

37N03M031002 

3 


MX 

LAPORTE 

37N03WQ3i003 

3 


BR 

LAPORTE 

37N03H03100A 

4 


RR 

LAPOR TE 

37N03W03H01 

3 


FH 

LAPORTE 

37N03H031102 

2 


FH 

LAPORTE 

37N03W031103 

3 


WT 

LAPORTE 

37N03W03110<i 

3 


WT 

LAPORTE 

37N03M031201 

1 


MX 

LAPORTE 

37N03H031202 

2 


FH 

LAPORTE 

37N03H031203 

3 


WT 

LAPOR TE 

3 7N03W03120A 

4 


BR 

LAPORTE 

37N03W031301 

4 


MX 

LAPQRTE 

37N03W03T302 

2 


FH 

LAPQRTE 

37N03W031303 

2 


BAA 

LAPORTE 

37N03W031304 

3 


RR 

LAPQRTE 

37N03W031401 

3 


MX 

LAPORTE 

37N03W031402 

3 


FH 


The procedure use to encode the raw data was described by Yahner, 
in “Using the Soil Survey for Land Assessment. A Computer Method". 
Minor modifications were made to Yahner' s source deck, in order to 
encode and manipulate the data necessary for the CZM inventory. 

Figure 6 shows an example of the output map and summary tables show- 
ing the acreage by potential productivity and potential erodibility. 

An addition was also made to Yahner 's programs that allowed us 
to select sample points for the field inventory. A discription, flow 
chard, program listing and tape file examples for program NEWSAMP 
appears in Appendix I. The output of NEWSAMP, fig 7, is a table show- 
ing the random grid locations for points in sections within a township. 
Points are selected with a probility proportional to the productive 
acreage. A summary at the top of the table shows the number of points 
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selected by productivity class. With this information, a topographic 
map, and a 1:24,000 grid the field crew could locate specific plot 
centers. 

Normal IDNR forest inventory procedures were followed for the 
field data collection. Five variable radius plots were measured for 
each sample point. The plots were arranged as a cluster sample, figure 
8. Tally trees were selected by a pps (probability proportional to 
size) method with a 10 factor basal area angle prism. Information re- 
corded on field sheets, figure 9, included: species, diameter at breast 

heighth (dbh), merchantable height (as number of logs), cull (percent 
of defect) and log grade. Ten year growth data for one dominant or 
co-dominant tree per points was measured. Also site index was measured 
for one dominant tree per sample plot. These data were sent to the IDNR 
offices in Indianapolis for card punching and analysis. 

The IDNR forest data processing service computer programs reduced 
the raw data to: 

• Stand and Stock tables, figure 10 

• Table of board foot volume by species 
and log grade, figure 11 

• Pulpwood production, figure 11 

• Statistical sample surmiary, figure 12 

• Ten year predicted stock and stand tables, 
figure 12 and 

• Average board foot growth per acre, figure 12. 

These results were produced for each sample point. Results for 

sample points were combined by Landsat classes to yield summary data 
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Figure 10. Example of Stand and Stock table output from field inventory data 
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Figure 11. Table of board foot volume by species and log grade, and summary 
pulpwood production table. 
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by spectral class. These data, supplemented by results from the 
previous years inventory form the basis for the recommendation of 
the study. 

RESULTS 

Evaluation of the inventory results indicate that forest lands 
in the Coastal Zone Management area: 

■ Predominantly belong to the Oak-hickory association 

• Are producing less average board foot volume per acre 
than their potential 

o Have been subjected to less than desirable management, 
and as a result 

• Have a majority of their growing fiber volume to low 
quality logs. 

• Are frequently growing on soils which are susceptible 
to erosion and, therefore, must be carefully managed. 

Table 1 sunmarizes the composition of the four forest classes 
which occur in the CZM. Typically the coastal zone's forests belong 
to the broad classification of forest ecosystem known as the central 
hardwoods. This association is composed of a mixture of hardwood 
species, which are predominated by the species oak. Rarely, do homo- 
genous stands of timber exist in this region. Usually, hardwood species 
intermix forming a large continua rather than a small uniform type. 
Certainly, the forest resources in the CZM reflect this intergration, 
as shown in figure 13. These diagrams represent the distribution, by 
volume Csawtimber trees), for the four classes of forest. 



Table 1. Description of Forest Classes in the CZM 


Class 


Symbol 


Sawtimber 


Pulpwood Bd.ft. growth/ 

ac»/year 


Oak-hickory- mapl e 

OHM 

3851 bd.ft.*/ac.vol. 
50 trees/ac.' 

54 sq.ft, basal area 
61% Oaks by vol. 

5 tons/ac.vol. 

79 trees/ac. 

22 sq.ft, basal area 

243 


2 

Oak-hickory- lowlands 

OH-l 

6019 bd.ft./ac.vol. 
55 trees/ac. 

70 sq.ft, basal area 
68% Oaks by vol. 

4 tons/ac.vol. 

64 trees/ac. 

16 sq.ft, basal area 

363 

CO 

CO 

2 

Oak-Hickory-Uplands ^ ^ 

4 . g£ 

Transition type 

OH 

5292 bd.ft./ac.vol. 
55 trees/ac. 

67 sw.ft, basal area 
78% Oaks by vol. 

3 tons/ac. vol. 

59 trees/ac 
15 sq.ft, basal area 

258 


T 

3823 bd.ft./ac.vol, 
43 trees/ac. 

51 sq.ft, basal area 
75% Oaks by vol. 

5 tons/ac.vol. 

72 trees/ac. 

20 sq.ft, basal area 

129 



* bd. ft. vol. by Doyle Scale. 
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The four types (fig. 13) belong predominantely to the oak-hickory 
association. The major differences occur in type one and four. The 
first type contains the greatest proportion of maple volume per acre. 

The fourth type has no maple volume but contains more red oak and 
cottonwood, so therefore, is considered a pioneering or transition 
type. Types two and three are more typical of the oak-hickory assoc- 
iation, each containing approximately equal amounts of the white and 
black oak groups and the walnut groups. The major difference between 
these types is that type two will more likely be found on moister lower 
sites than type three. 

For each type, the average annual productivity, measured by board 
foot growth per year, was calculated from the inventory data. These 
results are shown in the last column of Table 1. These productivity 
figures were compared to the potential board foot/acre productivity 
published by the Soil Conservation Service. 

Figure 14 shows the distribution of forest productivity classes 
in the CZM mapped according to potential board foot growth per acre 
per year. This data were summarized from Soil Conservation Service County 
Soils maps for Lake, Porter and La Porte counties. The forest acreage 
ranked by potential productivity from soil survey data is shown in Table 
2. A summary of Table 2 indicates; 

• 45% or 25,758 acres of forest land in the CZM is 
composed of pioneering species occupying sites 
capable of producing more volume and better quality 
timber than at present. 

• 39% or 23,191 acres of forest land in the CZM consist 



0 






Table 2. Distribution of Forest acreage by class for potential productivity 

GROUPS BASED ON ACCEPTED SOIL CONSERVATION SERVICE WOODLAND SUITABILITY 
RANKINGS, 





Potential Productivity 


cnorcT Tvoc 


High 

Medium-High 

Medium 

Low 


rUKto 1 Trt 








bd.ft ./ac/yr 

300- 

220- 

155- 

115- 



growth 

(Doyle Scale) 

375 

300 

260 

155 





(Acres) 



1 

Oak-hickory- pi e 

243 

1233 

2030 

1015 

145 

2 

Oak-hickory-lowiands 

363 

1668 

3626 

2321 

290 

3 

Oak-hickory-uplands 

258 

3843 

5656 

6671 

2103 

4 

Transition hardwood 

p 

129 

2683 

4714 

11530 

3553 
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of white, b.lack and red oaks, of small sawtimber 
dimensions growing on soils that have optimum 
site conditions for these species, and therefore 
are producing within acceptable productivity limits. 

• 15% or 8,919 acres of forest land in the CZM is com- 
posed of red oak, cottonwood and willow, which occupy 
sites that pose extreme management limitations, but 
are producing above average growth. 

Productivity, as used in Table 2, is a measure of a stands cap- 
ability to utilize the site (soil resource). The annual board foot 
growth does not suggest anything about stand quality or stocking. In 
order to determine occupancy (stocking) the average stocking per acre 
of Acceptable Growing Stock Trees (AGS) was calculated. These results 
are presented in Table 3. AGS trees. Table 4, are those species for 
which a timber market exists. Acceptable stocking for well managed oak- 
hickory stands is 60 to 100% in AGS trees. As Table 3 indicates, average 
stocking for forest types in the CZM is less than desirable. 

An evaluation of product quality further supports the contention 
that timber stands in the CZM are poorly managed. Table 5 shows the 
distribution of board foot volume by Purdue log grade. Greater than 
half the stand fiber volume is the study area is in poor quality logs. 

The high proportion of poor log grades indicates suboptimal stand 
management. This assumption is supported by comments from the field 
crews, which indicated that a large portion of the sample plots had 
been grazed, burned or high graded in the past. 

The data base created from the soil survey maps was manipulated 



TABLE 3. Growing Stock Distribution, in basal area and percent stocking. 



Acceptable GS 

Unacceptable GS 


Species 

Class 

Saw Timber 

Pul pwood 

Saw Timber 

Pulpwood 

Total 

OHM 






BA 

44.1 

16.8 

10.2 

4.9 

76.6 

% 

58 

22 

14 

6 


OHu 

• 





BA 

61.5 

12.6 

8:6 

3.5 

86.6 

% 

71 

15 

10 

4 


OH 






BA 

59.3 

11 .4 

7.7 

3.7 

82.6 

% 

72 

14 

9 

5 


T 






BA 

47.9 

14.0 

3.2 

6.2 

73.3 

% 

65 

19 

7 

9 



Acceptable Growing Stock are all merchantable tree species. Table 4. 

Unacceptable Growing Stock trees are those species for which presently no market exists. 
Saw timber are trees which d.b. h. is greater than 10 inches. 

Pulpwood trees are those between 6 and 10-inch d.b.h. 

BA - basal area is a measure of stem occupancy, expressed in square feet/acre. 
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TABLE 4. Acceptable Growing Stock trees for types in the CZM. 


White Ash 

Sugar Maple 

Basswood 

Soft Maple 

Beech 

White Oak 

Cottonwood * 

Red Oak 

Cherry 

Black Oak 

White Elm 

Tulip Poplar 

Shagbark Hickory 

Sycamore 
Black Walnut 


TABLE 5. Distribution of board foot volume by log grade for forest 
types in the CZM. 


Type 


Purdue Log Grade 

Prime 1 2 ; 

(percent of board foot vol/ac) 


Oak Hickory Maple 

8 

9 

33 

50 

Oak Hickory lowlands 

13 

13 

31 

43 

Oak Hickory uplands 

10 

11 

30 

49 


Transition 


5 


6 


23 


66 
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to give generalized maps for; 1) potential board foot/year produc- 
tivity (fig. 14), and 2) soil capability (fig. 15). The soil cap- 
ability map, figure 15, indicates generalized zones where soils possess 
similar characteristics, related to their susceptibility to erosion. 

Zone I is the least susceptible. Zone III the most susceptible and, 
therefore, will require the most caution when planning for any human 
activity on these soils. 

Table 6 gives the distribution of forest area by soil capability 
zones. 

RECOMMENDATIONS 

The recommendations in this section are based on the preceeding 
results and accepted management practices for the oak-hickory cover 
type, guidelines are given for water, wildlife, recreation and timber 
uses of the forest resource. 

Water 

Since the majority of forest land in the coastal zone occurs on 
erodible soils, these areas should be a carefully managed to prevent 
excessive erosion. Human activities, including the development of hous- 
ing sites, should be regulated in order to minimize soil loss and 
stream sedmentation. 

Intensive silvicultural activities should be limited in Zone III 
(figure 15) forest lands. Only those activities that perpetuate good 
watershed management should be utilized. The management of these 
forest lands for a commercial timber crop should be carefully monitored. 
In no instance should indiscriminate timber clearing be condoned. 


/ 
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Figure 14;, 

Potential Forest Productivity 


.+++4 



High 300-375 bdft/ac/yr 
Medium High 220-300 bdft/ac/yr 
Medium 155-260 bdft/ac/yr 


Low 115-155 bdft/ac/yr 
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Wildlife 

Wildlife management in the oak-hickory forests is compatible 
with other management practices. These forests provide habitat, for 
numerous species. The enhancement of wildlife habitat usually occurs 
in conjunction with selective cutting or timber stand improvement 
practices. These practices produce small openings in the forest stand 
that encourage wildlife. 

If wildlife is a highly desired element in a management scheme, 
specific activities should be pursued. Ideally, a manager should strive 
for a stand having a good distribution of size classes from saplings to 
mature sawtimber. Currently, forests of the coastal zone lack the ideal 
with regard to diameter distribution. However, through controlled thin- 
nings and timber stand improvement methods stand structures could be modi- 
fied. In order to make these activities economically feasible, these 
practices should be concentrated to forest areas growing on highly pro- 
ductive soils (figure 14). 

Recreation 

The coastal zone management area is located within easy reach of 
millions of people living around the fringe of Lake Michigan. Recreational 
opportunities are currently at a premium for such a large populace. The 
coastal zone forests along the lake and major water courses offer the 
potential for recreational development, 

Figure 16 delineates areas where recreation is best suited for 
development. However, care should be taken when developing these areas 
to prevent misuse of the soils. The forests along the coast grow on 
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TABLE 6. Distribution of forest land by 

Zone Occupancy 

(percent forest) 

I 40.5 

c 

II 4.0 


III 


55. 5 


CAPABILITY IN THE CZM . 


Description 


No limitations due to excessive 
soil erosion or slope. Usually 
wel 1 dra 1 ned soil s . 

' tn 

o 

Limitation because of slope or wet- 
ness and/or erosion hazards are more 
serious, and require care in manage- 
ment . 


These soils have severe limitations, 
often due to excessive slopes and 
generally unstable soil profiles. 
These areas should have use restrict- 
ed to activities which do not promote 
removal of the vegetation. 
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unstable soils prone to excessive wind and soil erosion. In these 
areas summer residences and large camping complexes should be kept 
to a minimum. The area is best suited to smaller recreational sites 
of a dispersed nature. 

The area along Deep River should be managed for maximum protection 
of the stream course. Heavy development should not be permitted along 
the stream. Walking trails or camping/picnic facilities should be 
located with great care. 

Timber 

The potential for developing the timber resources of the coastal 
zone is best for areas growing on highly productive and least erodible 
site. These sites should be managed for forest production (fig. 16). 
Currently these areas are not optimally managed. Practices which 
would increase productivity in these stands should be encouraged. 

Basically, silvicultural practices, primarily thinnings with some 
selective cuttings should be undertaken on these sites. These activ- 
ities should be pursued with the objective of perpetuating oak regen- 
eration and increasing stocking of the acceptable growing stock trees. 

If intensive management activities were pursued over a 40-or 50-year 
planning horizon an estimated 56 percent increase in timber value (using 
current timber price information) could be anticipated. 

Currently, limited opportunities exist in the coastal zone to fully 
exploit the timber potential. Figure 17 shows the number of sawmills in 
the counties adjoining the study area. This lack of primary wood using 
industry in the CZM creates difficulty in marketing timber for top prices. 
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Figure 17. Number of sawmills in counties adjoining the CZI'^. 
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CONCLUSIONS 

This report culminates one phase of a cooperative program 'between 
LARS and the IDNR- Division of Forestry. The results of this two year 
effort, documented here and in previous semi-annual reports, can be 
summarized as follows: 

• Computer-aided analysis of Landsat digital 
data can be effectively used to help 
characterize the central hardwood forest 
resource. 

• Landsat classification provide timely 
information about the forest resource base 
which is not available elsewhere 

• Spectral differencies of forest classes 
classified on Landsat data can be used as a 
basis to allocate field sample plots, and 
therefore, reduce costs associated with 
collecting the field data 

• Landsat classifications in combination with 
ancillary information, such as soils maps 
provide a powerful tool to help the IDNR set 
priorities for forest management opportunities. 

Empodied in the recommendations is the knowledge that use of remote 
sensing (i.e.: computer-assisted Landsat classifications) can benefit 
man. These benefits accure to the populace indirectly through legisla- 
tive actions that protect and improve Indiana's forest resources. The 
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CZM study results and recommendations have identified the potential 
for various management alternatives. These recommendations will be 
translated into land use policies. Empowered by federal mandate, 
the State of Indiana will enforce these policies through zoning regu- 
lations designed to protect the coastal zone management area resources. 

Landsat inputs, although minimal, form a critical base upon which 
the forest resource policies will develop- Without this remote sensing 
input the quantative data upon which to base this studies recommendations 
would not have been available. Remote Sensing has definately played an 
important role in shaping Coastal Zone Management policies for the forest 
resources in northern Indiana. 
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Indiana Coastal Zone Management Project 
Summarization and Sampling Program Documentation 


PURPOSE OF THE PROGRAM 


The name of this program is "SAMPLE". Its purpose is to 
summarize and sample from data collected on a township basis. 
Given the soil type of a cell with forest cover, the program 
assigns that cell a value for the soil's productivity, and a 
value for its capability. It then counts the total number of 
acres per productivity class and per capability unit on a sec- 
tion basis and a township basis. The n\amber of cells per pro- 
ductivity class which need to be field-checked are calculated, 
and then are randomly selected from the township. See the 
following pages for more details. 


THE EXEC ROUTINE 


This program is written for a IBM 370 computer in fortran 
IV, G level, in conjunction with a CMS 360 exec routine. This 
exec routine has a file name of SAMPLE and a file type of EXEC. 
This means that whenever "SAMPLE" is typed on CMS mode, state- 
ments are defined and commands are carried out which are found 
in the file SAMPLE EXEC. 

Contents of the exec file: 


GLOBAL TXTLIB SSPLIB SKSLIB 


FILEDEF 
FILEDEP 
FILEDEF 
FILEDEF 
FILEDEF 
FILEDEF 
FILEDEF 
GETDISK TEMP 
LOAD SAMPLE 
START SAMPLE 
RELEASE 192 T 


DSK LAKE SOIL 
DSK PORTER SOIL 
DSK LAPORTE SOIL 
DSK TWNSHP &1 
DSK SEED &2 
PRT 

8 DSK-Tl LRECL 10240 ETENT 
5CYL CLEAR 


(DET 


36 

ORIGINAL PAGE IS 
OF POOR QUALiry 


• The statement which begins the exec routine enables the 
program to have access to the scientific sub-program library, 
which IS used to generate random numbers in this case. Access 
to the system library is also obtained. 

Next follow seven statements which define program files: 

Program file 1 is defined as the file which is on disc 
with the file name and file type of LAKE SOIL. Program file 
2 is defined as the disc file PORTER SOIL, and program file 3 
is defined as the disc file LAPORTE SOIL. 
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For the program files 4 and 5^ another note of explanation 
needs to be made. In exec routines, an ampresand (&) may be 
used which will be exchanged for a value later on. This ex- 
change is made when you are in CMS, and "SAMPLE", followed by 
a space and a number, is typed on the terminal. For example, 
with the exec statement: 

FILEDEF 4 DSK TWNSHP &1 

and when "SAMPLE 369" is typed in the CMS mode, the exec rou- 
tine will be performed and the program file 4 will be defined 
as the disc file TWNSHP 369. ‘ 

More than one value may be exchanged. The ”1" placed 
after the ampersand in the FILEDEF 4 statement indicates that, 
in CMS, the first value typed after "SAMPLE" is typed will be 
exchanged for the &1, and the second value typed after "SAMPLE" 
will be exchanged for the & with the 2 placed after it, as in 
the FILEDEF 5 statement. The first and second values are se- 
parated by a space. When "SAMPLE 369 5" is typed in CMS, the 
exec routine will run with the program file 5 defined as the 
disc file SEED 5. 

Program file 6 is defined as the printer, which means that 
when the program says to write things on file 6, the information 
will be printed as output. 

Program file 8 is defined as a temporary disc with a speci- 
fied format. This file also must be defined within the program 
with a fortran statement. 

The exec routine also gets a temporary disc of five cylin- 
ders, clears all previous information and logs xt into your 
system. 

The next commands load for fortran program with the filename 
of SAMPLE, and begin its execution. Once the execution of the 
program is completed, device 192, which is the terapory disc, is 
detached from your system. 

Therefore, to execute the program while CMS, simply type in 
SAMPLE followed by the TWNSHP file type you want, and with the 
SEED file type you want. 


FORMAT OF INPUT 


1) . TWNSHP 

There are five sources of input to this program. All are 
located on disc. The major data which the program works on is 
contained in the file with the name TWNSHP and of the type which 
is a series of three integers, depending on the particular TWNSHP 
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for which you want the output. For example, a file called 
TWNSHP 369 means that the program will work on the data from 
township 36 N, range 9W. The other TWNSHP files follow a 
similar pattern. 

The TWNSHP file contains the following information in the 
following format, all on one line. Each coordinate row and col- 
umn combination which has data is on a separate line. The 
number in parenthesis is the number of spaces which the data 
may jbake, or the number of spaces between pieces of data. 

County name (8), Township (3), Range (3), Section (2), Coordin- 
ate row (2), Coordinate column (2), Forest spectral class (2), 
(6), An end of section marker (3), Soil Symbol and slope (3), 
Erosion value of soil (1) . 

The end of section marker is either the letters ALL or HAF. 

When data for an entire section is available, use ALL. When 
data for less than an entire section is available, use HAF. 

See page 5 for an example of the data in a TWNSHP file. 


How data for a TWNSHP file was collected : 

The data for a township was encoded on a section basis. 
Coded LANDSAT data for the section was used . Each section was 
divided into 256 2%-acre cells. A particular cell's location 
was defined by a coordinate system consisting of 16 rows and 
16 columns, which also totals 256 cells. Information for the 
section was collected on each cell which had a forest spectral 
class in it, which included the spectral class, and the soil 
type found at that location. The soil type was found by using 
a soil map. 


2. Soil file 

For reasons of convenience, much of the information having 
to do with the soil type was not put into the TWNSHP file, but 
was put into different soil files, depending on the county in 
which the soil was located. The program matches this informa- 
tion with the soil types in TWNSHP. 

The soil files have the name of either LAKE SOIL, PORTER 
SOIL, or LAPORTE SOIL. The following information is contained 
in these files, in the following format. The number in parenthesis 
is the number of spaces which the data may take on a line, or the 
number of spaces between pieces of data. 

Soil symbol and slope (3) , soil name (6) , soil texture (6) , 
erosion number (1) , soil productivity class (1) , soil capability 
unit (1) . 

See pages 6-9 for a copy of the information which is in the soil 

ORIGINAL PAGE IS 
OP POOR QUAUTYi 
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3) . SEED 

The last file which is used as input is the file named 
SEED. The sub- program which calculates the random numbers 
calls for a seed with which to begin its calculations. This 
number should be a prime number. To avoid using the same 
random numbers for each TWNSHP file which is run through the 
program/ a different seed is used for each TWNSHP file that 
is run. These each are simply contained within files with the 
name SEED and of type designated by a number. The file type 
simply differentiates one seed's file from another. Each 
file contains one seed. 


HOW THE PROGRAM SAT4PLES 


The program gives the location of a particular number of 
235 -acre cells which should be field-checked to find out what 
is actually growing in that cell . 

When the data ^f or each section is read and complied, it is 
stored on a temporary disc. The entire section's data takes up 
one line on the disc. Since there may be up to 36 sections in- 
cluded in a township, there may be up to 36 lines on the tempor- 
ary disc. The section's data within the line where it is stor- 
ed is made up of the following information for each 23s-acre cell: 

The row coordinate of the cell 
The column coordinate of the cell 
The forest spectral class form the Landsat data 
Whether data for the entire section was available 
or not, 

{i.e. only entered for the last 2H acre cell of 
the section. ) 

The soil symbol and slope 

The soil ' s erosion value 

The soil name 

The soil texture 

The soil's productivity class 

The soil's capability unit 

The section number. 

I 

Next, the number of 235 -acre sample cells required per pro- 
ductivity class is calculated. This is done by taking 1%' of 
the number of 235 -acre cells which are in each productivity class. 
This number is rounded off to the nearest integer. 

400 groups of random numbers are then called. The first 
two random numbers of the group are between 1 and 16, inclusive, 
and the third is between 1 and 36, inclusive. When a group of 
3 random numbers is called, the third number points to a row on 
the temporary disc which contains the township information. That 
is, it points to a particular section's data. The first two 
random numbers then pick out a row and column coordinate location 
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of a 2%-acre cell in the section. If the soil productivity 
class of that cell is one from which a sample is required, 
that cell location, etc. is printed. If it is not, another 
group of 3 random numbers is tried. This is repeated until all 
the samples which are required have been found. Since there 
are many cells in a section with no forest spectral class and 
hence no productivity class, this procedure for finding the 
sample cells required may be repeated frequently. Since there 
are only 400 groups of random numbers, this procedure may not 
be repeated more than 400 times, however. Hopefully by then 
all the samples will have been found. 


PROGRAM OUTPUT 

The following is printed in the output: 

For each section - County, Township, Range, Section number 

An indication for whether all the data for 
section was available or not. (If all 
not available, an asterisk is beside the 
section no . ) 

No. acres per productivity class 
No. acres per capability unit 
, Total no. of acres in all productivity classes 
Total no. of acres in all capabilxties units 

For the township - No. acres in each prod, class for all 

sections . 

Total no. acres in all prod, classes for all 
sections . 

No. acres in each cap. unit for all sections. 
Total no. acres in all cap. units for all 
sections . 

The no. of sections for which all the data was available minus 
the no. of sections for which not all the data was available 
(called a half-section) . 

Total no. of sections for which data was taken. 

The no. of samples required per productivity class. 

Total no. of samples required. 

The seed which was used for the random niimber calculations. 

For the sample cells - County 

Township 

Range 

Section 

Row Coordinate 
Column Coordinate 
Soil symbol and slope 
Soil erosion value 
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Any "misspellings" of soil symbols in the TWNSHP file, with 
the section, row, and column of the misspelling. 

The nuinber of soil names (NUMB) in a county's soil file if 
there are more than 99. 







ORIGINAL PAGE 
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Wrj.te the sums 
in a table 



Write table headings headings 
for location of samples 


Write no. sections 
read 



Randomly pick one of the 
sections written on the diski 
with the third random no 
in the triplet. 



/ A prod, class at a \ 
random row+col - e prod, 
class that must be 
sampled from j/ 

7 X 


Decrease no. of samples re- 
quired for that prod, 
class 



"Write. 

County ROW 

TNSP COL 

RMG soil symbol 

Section soil slope 

erosion 
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After 400 X around 
the loop. 


STOP 
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'ILE. . . NEWSAMP fortran PI 

*»e«eSUMMAHI7ATlOAI AND SAMOLljgG PR0'3R AM**^'^* 

WRITTEN BY MARK BEAN 
JUNE 28, 1977 

INDIANA COASTAL 20NL MANAGEMENT PROJECF 

GLOBAL TXTLIB SSPLIB SYSLIB 
FILEDEF 1 DSK LAKE SOIL 
FILEDEF 2 OSK PORTER SOIL 
FILEDEF 3 OSK LAPORTE SOIL 
FILEDEF 4 DSK TWNSHP Nl 
FILEDEF 5 DSK SEED a.2 
FILEDEF 5 PRT 

FILEDEF 8 OSK-Tl LRECL li)240 BLKSIZE 10240 XTENT 36 

GETDI5K TEMP 5CYL CLEAR 

LOAD SAMPLE 

START SAMPLE 

RELEASE 192 T (DET 

ENTER 'ALL' IN COLUMNS 30-32 WHEN DATA FOR AN ENTIRE SECTION 
HAS BEEN ENCODED. 

enter 'HAF* in COLUMNS 30-32 WHEN HALF OF THE DATA FOR A SECTION 
HAS BEEN ENCODED, AND NO MURE DATA FOh THIS SECTION 
IS AVAILABLE. 

Integer soiluoo,6> ,sectn(16.i6,io) ,nsa’ip(6) ,sel,uow,col,all,haf, 

1 BLANK, HOLD, PORT 

DIMENSION ACC(36,9) ,ACR(36,9) ,UNCC36,9) ,UNP(36,9) ,SAMP((j) , 

1 NX(400) ,NY(400) ,NZ(400) ,TUwP(9) ,TACP(9) ,TACC(9) 

DATA ALL, BLANK, HaF , LAKE , PORT . LAPO /'ALL ' • , »HaM, 

1 'LAKE' , 'PORT* , *LAP0»/ 

REAL os CNTY 

integer "4 key 

EQUIVALENCE (CNTY, KEY) 

DEFINE file 8<36,2560,‘J,IPTR) 

IPTR=1 

HNO=0. 

TNO=0. 

NN1 = 1 
NN2 = 1 
NN3=I 
NO=l 

INITIALIZE AMO read a card 

5 CONTINUE 
DO 20 1=1,16 
DO 25 J=l,16 
DO 30 K=l,10 
SECTNd, J,K) =0 
30 CONTINUE 
25 CONTINUE 

20 CONTINUE 

DO ?1 1=1,100 
DO 22 J=l,6 
50IL(I,J) =0 
22 CONTINUE 

21 CONTINUE 

12 HOLD = BLAN^ 

READ (4,300, END=140) CNTY, TNSRtRNG, SEC, ROW, CO l,SECTN(HO„, COL, 1 ) , 

1 HOLO,(SECTN(ROw,COL,N) ,N=3,A) 

300 FORMAT(A8,2A3,3I2,6a,I3,A3,6X,A3,I1) 

SECTN(R0W,C0L,2) = MOLD 
IF(KEY.EQ.LAKE)G0 TO 6 
IF(KEY.Ea,PORT)GO TO 7 
IF (KEY.EQ.LAPO)GO TO « 

READ THE NO. OF SOIL TyPE'=. COUNTY NAME, ANU SOIL DATA 


OF POOR QUAUIYI 
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NEWSAMP FORTRAN PI 


C PEP00790 

6 REWIND 1 P|£92®9S 

IF(SECTN(R0W»C0L>2) .NE. ALL. AND. NNl.GE. 3)00 TO 13 5IES9S4S 

IF <SECTN(ROW»COLf2) .NE.HAF.AND.NN1.GE.2)G0 TO 13 PEP00B20 

read(1.ioo)numb 

100 FORMAT (13) 

IF(NUMB.GT,99)G0 TO 980 PEP00850 

READ(1.200) t (SOIL(ItJ) *3=1.6) .1 = 1, NUMB) PEP00B60 

200 F0RHAT(A3,2A6,3I1) '’gPOOglS 

NN1=NN1+1 PEP00880 

GO TO 13 PEP00890 

7 REWIND 2 PEP00900 

IF(SECTN(R0W.C0L.2) .NE.ALL.AND.NN2.GE.2)G0 TO 13 PEP00910 

IF(SECTN(R0W.C0L,2) .NE.HAF.AND.NN2.GE.2)G0 TO 13 PEP00920 

READ(2.100)NUMB PEP00930 

IF(NUMB.GT.99)G0 TO 980 PEP009AO 

READ (2,200) ((SOIL(I.J) , J=1 , 6) , 1=1 , NUMB ) 

NN2=NN2-.l 

GO TO 13 PEP00970 

8 REWIND 3 PEP00980 

IF(SECTN(R0W,C0L.2) .NE.ALL.AND.NN3.GE.2)G0 TO 13 

iF(SECTN(R0W,C0L,2) .NE.HAF.AND.NN3.6E.2)G0 TO 13 PEPOIOOO 

READ(3,100)NUMB ElSSlSiS 

IF(NUHB.GT.99)G0 TO 980 PEP01020 

READ(3,200) ( (SOIL(I.J) , J=1 ,6) , 1=1 .NUMB) 

NN3=NN3+1 PEPOIOAO 

GO TO 13 £159^919 

C PEP01060 

c pEPOioao 

C SET THE SOIL DATA INTO THE SECTION ARRAY PEP01090 

c pEpOiioo 


13 IF(ROW.EQ.O)GO TO 80 
INUMB=NUMB+1 
DO 50 M=1.INUMB 

IF(SECTN(ROW,COL,3) .NE.SOIL(M.I) )GO TO 60 
SECTN(ROW,COL,5)=SOIL(H,2) 

SECTN ( ROW , COL . 6 ) =SOIL ( M , 3) 
SECTN(R0W.C0L.8)=S0IL (M.5) 
SECTN(R0W,C0L,9)=S0IL(M,6) 

GO TO 52 

60 IF(M.EQ.1NUMB)G0 to 55 

51 CONTINUE 
50 CONTINUE 

52 INUMB=INUMB-I 
SECTNtROW.COL.lO) = SEC 
IF(SECTN(R0W,C0L,2) .EQ.ALL) GO TO 80 
IF(SECTN(R0W,C0L,2) .EQ.HAF) GO TO 80 
GO TO 12 


PEP01120 

PEP01130 

PEP01140 

PEP01150 

PEP01160 

PEP01170 

PEP01180 

PEP01190 

PEP01200 

PEP01210 

PEP01220 

PEP01230 

PEP01240 

PEP01250 

PEP01260 

PEP01270 

PEP01280 

PEP01290 


C PEP01300 

ro«-tt**-tt**tt*-it**«-*»it***-iK*«e**«o««»o*-t**«*-«-**'tt«««*<*«<.«-^****»*'<f*«*<H>*******»»PEP0131 0 
C PEP01320 

C INITIALIZE AND ACCUMULATE PEP0j330 

C (SECTION BASIS) PEP01340 

C PEP01350 

r»0» «•*«»<*****> **<*<*» ****<>«**<><*'►* ***PEP 01 36 0 

c PEP01370 

80 CONTINUE PEP01360 

NN1=1 PEP01390 

NN2=1 PEP01400 

NN3=1 PEP01410 

IF(SECTN(R0W,C0L,2) .EO.HAF) HnO = HNO + 0.5 PEP01420* 

DO 90 J=l,9 PEP01430 

ACC(NO,J)=0. PEP01440 

ACP(NO,J)=0. PEP01450 

UNC(NO,J)=0. PEP01460 

UNP(N0.J)=0. PEP01470 

90 CONTINUE PEP01480 

DO 110 1=1.16 PEP01490 

DO 115 J=l,16 PEPOISOO 

DO 120 K=1.9 „ PEP01510 

IF(SECTN(I, J.S) .EQ.K)UNP{NO,K)=UNP(NO,K) +1. PEP01520 

IF (SECTN ( I, J, 9) .EQ.K)UNC(NO.K)=UNC(NO,K) +1. PEP01530 

120 CONTINUE , PEP01540 

115 CONTINUE / PEP01550 

no CONTINUE PEP01560 
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FILE. . . NEWSAMP FORTRAN PI 

C PEP01570 

PEPOlS'^O 

CALCULATE ACRE*;, TOTAL, AND rfHlTE PEP01600 

PEPO i 61 0 

PEP01630 

TOTC=0. 

TOTP=0. PEP016b0 

DO 130 J=l,9 PEP01660 

ACP(NO,J)=UNP(NO,J)«2.‘5 PEP01670 

ACC(NO,J1=UNC{NO,J)*2.S PEP01680 

TOTP=sACP(NO,J) +TOTP PEP01690 

TOTC=ACC(NO,J)+TOTC PEP01700 

130 CONTINUE PEPOl^O 

IF (NO.EO.DfjO TO 106 ' PEP01720 

IF{5ECTN (R0vil,C0L,2) .EQ.HAF) (3 o TO 107 PEP01730 

GO TO 103 PEPOITA^O 

_ _ PEP017b0 

PEP01760 
PEP01770 
PEP0I780 


5 uRITE(6,6A0)CNTy,TNSP,RNG PEP017b0 
) F0RHATtlHl//l5X, 'COUNTY AS . 3 X TOWNSHIP A 3 . bX RANGE ',A3, PEP01760 
1////17X, 'ACRES PER PRODUCTIVITY CLASS', 24X, 'ACRES PER PEP01770 
2»CAPABILITY UNIT ' ,/lX, ' SECTION' , AX 1 * , 6 X, »2' , 6 X, ' 3' , PEP0I7S0 
36X, ' A' , 6 X, *5* , 6 X, ' 6 ' ,9X, ' 1' , 6 X, *2* , 6 X, '3' , 6 X, ' A' , 6 X, *5' , 6 X, ' 6 * ,bX» PEPO 1790 


A*7'|6X,'8',6X,'9') 

107 WRITE <^?lo3) SEC, ( ACP ( NO , I ) , I =1 , 6) , ( ACC (NO, J) , J=i , 9) , TOTP 


PEPoiaoo 

PEPojaio 

PEP01S20 


601 FORMAT{/3X,I2, "" ,AX,Fb.l,2X,FS.l,?X,F5.1,2X,Fb.l,2X,F5,l,2X,F5.1 PEP01S30 

1,AX,9(F5. 1,2X) ,AX,Fb.l) PEP018A0 

GO TO 109 PEPOlSbO 

108 WRITE (6,602) SEC, ( ACP ( NO , I) , I = 1 , 6 ) , ( ACC ( NO , J) , J=1 , 9) PEPOlSfaO 

602 F0RMAT(/3X»I2,SX,F5.1,2X,F5.1,2X,F5.1,2X,F5.1,<;X,F5.1,2X,F5.1,AX, PEPO 1870 
19(F5,1,2X) ,AX,F5.1) 

109 WRITE(8'IPTR)SECTN PEP01890 

NO=NO+l PEP01900 

GO TO 5 PEP01910 

PEP01920 
PEP 0 1 93 0 
PEP01940 

INITIALIZE, ACCUMULATE, AND WHITE PEP019b0 

(TOWNSHIP BASIS) PEP0I9b0 

PEP01970 

»<n, »«*PEP 01 98 0 

PEP01990 

lAO NO=NO-l PEP02000 

DO 150 1=1,9 PEP02010 

TACP(I)=0. PEP02020 

TACC(I)=0. PEP020J0 

TUNP(I)=0. PEP020A0 

150 CONTINUE PEP020b0 

TTACP=0. PEP02O0O 

TTACc=o. 

DO 160 J=1.9 PEP02080 

DO 165 1=1, NO PEP02090 

TACP t J) =ACP (I ♦ J) +TAC° ( J) PEP0210 0 

TACC(J)=ACC(I ,J)+TACCfU) PEP02110 

TUNP(U)=dNP(I,J)-,TUNP(J) PEP02120 

165 CONTINUE PEP021J0 

160 CONTINUE PEP02U0 

DO 166 J=l,9 PEP021bO 

TTACP=TACP ( J) +TTACP PEP02160 

TTACC=TACC(J)+TTACC PEP02170 

166 CONTINUE PEP02180 

WRITE(6,6A1)CNTY,TNSP,PNG ^ PEP02190 

6A1 FORMAT (1H1//15X, 'COUNTY * , A8, 3X, ' TOMNSHIP A3 , 5X ,' RANGE *,A3) PEP02<^00 

WRITE(6,650) (K,TACP(K) ,K = 1.6) PEP02210 

650 F0RMAT(//////////20X, 'TOTAL ACRES PEk PRODUCTIVITY CLASS', /20X, PEP0222(; 

1 3A ( 1H=) //25X, ' PROD. ' , I AX , 'ACRES' . /?5a. 'CLASS • ,// 6 ( 27X , 1 1 , IbX, F6 . 1PEP022 jO 

PEP022‘fU 
PEP02250 
PEP02260 
PEP02270 


2 /) ) 

WR1TE(6,660)TTACP 
660 FORMAT(25X,'TOTAL»,12X,F7.1) 
WRITE (6,670) (K.TACCCK) ,K = 1,9) 


670 FORMAT (/////////21X, 'TOTAL ACRES PEP CAPABILITY UNIT • /2 1 X , 31 ( 1H= ) /PEP02280 
I /27X, 'CAP. ' ,12X, 'ACRES' ,/27X. 'UNIT* , //9 ( 28X , 1 1 , 1 3X , F6 . 1 /) ) PEP02290 

WRITE (6,680) TTACC PEP02J00 

680 F0RMAT(26X, 'TOTAL' ,10X,F7.1) PEP02310 

TNO = FLOAT (NO) - HNO PEP023c0 

WRITE(6,605)TNO,NO PEP02330 

605 FORMAT (/////5X, 'NO. OF SECTIONS MINUS HALF SECTIONS = * , FA. 1//5X , ' PEP023A0 


ORIGINAL PAGE IS 
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190 

690 


NEUSAMP FORTRAN PI 
TOTAL NO. OF SECTIONS = *»I2) 


NSAHbQ 

DO 190 Isl.6 
SAMPm=TUNP(l)*,Ol 
NSAMP<I)«>SAMP(I)*0.5 
NSAHoNSAHi’NSAHPa) 
CONTINUE 


FILE. . , 

PEP02350 

C PEP02360 

C PEP02380 

C CALCULATE NO. OF SAMPLES REQUIREDt AND PRINT PEP02390 

C PEP02400 

PEP02420 
PEP02430 
PEP02440 
PEP02450 
PEP02460 
PEP02470 
PEP024S0 
PEP02490 


WRITE<6|690) (I.NSAMPd) .1 = 1,6) ,NSAM 

. F0RMATtlHl//22X. 'NO. OF SAMPLE CELLS PER PRODUCTIVITY CLASSV/29X.PEP02500 

1 'PROD. •. iSX, • SAMPLES * //6 ( 31X, II »20X, 12/)/. 30X, • TOT AL*.16X, 13) PEP02510 

C PEP02520 

C PEP02540 

C CALL 400 RANDOM NOS. IN GROUPS OF 3 PEP02550 

C PEP02560 


READ <5,400)1X 
400 FORMAT (15) 

WRITE(6,405) IX 
405 FORMAT (/5X, I SEED = 
DO 210 1=1,400 
CALL RANDU(IX,IY,YF) 
IX=IY 

NX(I)=YF»16+1 

CALL RANDUdX.IY.YF) 

IX=lY 

NYd)=YF*16*l 
CALL RANDUdX.IY.YF) 
IX = IY 

NZd)=YF*NO+l 
210 CONTINUE 


• 15) 


PEP02580 

PEP02590 

PEP02600 

PEP02610 

PEP02620 

PEP02630 

PEP02640 

PEP02650 

PEP02660 

PEP02670 

PEP02680 

PEP02690 

PEP02700 

PEP02710 

PEP02720 

PEP02730 

PEP02740 


o P2P 02760 

C PEP02760 

C SAMPLE THE TOWNSHIP, DECREASE NO. OF SAMPLES REQUIRED. PEP02770 

C AND PRINT PEP02780 

C PEP02790 

C 

WRITE(6.689) 

DO 220 1=1,400 
IPTR=NZd) 


READ (8' IPTR)SECTN 

IF(SECTN(NXd) ,NYd) .8) . EO. 1 . AND .N5AMP (1 ) .GT . 0 ) GO 
IF(SECTN(NXd) ,NYd) .8) .EQ.2.AND.NSAMP (2) .GT.0)GO 
IF(SECTN(NXd) .NYCI) .8) .EQ.3. AND.NSAMP (3) .GT. 0) GO 
IF(SECTN(NX(I) ,NY(I) .8) .EQ.4. AND.NSAMP (4) ,6T.O)60 
IF(SECTN(NXd) ,NY(I) .8) .EQ. 5. AND.NSAMP (5) .GT.0)GO 
lF(SECTN(NXd) ,NY(I) .8) .EQ.6. AND.NSAMP (6) .GT.O)GO 
GO TO 219 

222 NSAMP(1)=NSAMP(1)-1 

WRiTE{6.691)CNTY,TNSP.RNG,SECTN{NXd) .NYd) ,10) ,N) 
l.SECTN(NXd) ,NYd) ,3) .SECTN(NXd) ,NYd) ,4) 

GO TO 219 

223 NSAMP(2)=NSAMP(2)-1 

WRITE(6,692)CNTY.TNSP,RNG,SECTN(NXd) ,NYd) .10) ,N) 
l.SECTNlNXd) ,NYd) ,3) .SECTN (NX d ) ,NY d ) . 4) 

GO TO 219 

224 NSAMP(3)=NSAMP(3)-1 

WRITE(6.693)CNTY.TNSP.RNG,SECTN(NXd) ,NY(I) ,10) ,N> 
1, SECTN (NX d) ,NYd) ,3) . SECTN (NX d ) , NY (I ) ,4 ) 

GO TO 219 

225 NSAMP(4)=NSAMP(4)-1 

WRITE(6,694)CNTY.TNSP,RNG,SECTN(NXd) ,NY(I) ,10) ,N> 
1, SECTN (NX (I) ,NY(I) .3) . SECTN (NX d) ,NY(I) ,4) 

GO TO 219 

226 NSAMP(5)=NSAMP(5)-1 

WRITE(6,695)CNTY.TNS'’,RNG,SECTN(NX d) ,NY ( I ) , 1 0 ) ,N) 
1. SECTN (NX (I) ,NY(I) .3) .SECTN (NX (I ) ,NYd) ,4) 

GO TO 219 


TO 

TO 

TO 

TO 

TO 

TO 



PEP02810 

PEP02820 

PEP02830 

PEP02840 

PEP02850 

222 

PEP02860 

223 

PEP02870 

224 

PEP02880 

225 

PEP02890 

226 

PEP02900 

227 

PEP02910 

PEP02920 

PEP02930 

.NY (I) 

PEP02940 

PEP02950 

PEP02960 

PEP02970 

.NYd) 

PEP02980 

PEP02990 

PEP03000 

PEP03010 

.NY (I) 

PEP03020 

PEP03030 

PEP03040 

PEP03050 

.NYd) 

PEP03060 

PEP03070 

PEP03060 

PEP03090 

,NYd) 

PEP03100 

PEP03110 

PEP03120 



III-5 


FILE. . . 


227 


NEWSAMP FORTRAN PI 


NSAMP(6)=NSAMP(6)-1 

WR1TE(6.696)CNTY.TNSP.RNG.SECTN(NX(I) ,NY ( I ) . 1 0 ) .NX ( I ) i NY 
1,SECTN(NX(I) .NY(I) rJ) . SECTN (NX ( I ) ,NY ( I ) .4) 

GO TO 219 

689 format (/////SIX. 'COORDINATE LOCATION OF SAMP1.ES' /SIX. 30 ( 
11 OX. 'COUNTY' *8X.' TOWNSHIP • .8X. 'RANGE' .8X. 'SECTION' ,8X.*R 
28X»'C0LUMN' »8X. 'PROD. CLASS '. 7X .' SOIL SYMB.'//10X. 


(II 


,H=)// 

IW' . 


691 

692 

693 

694 

695 

696 

219 

220 


3^14.(.1H-)///) 


. 0RMAT(11X.A8,7X.A3.11X.A3.11X.I2.12X.I2.11X»I2.15X.'1'.16X.A3.//) 
F0RMAT(11X.A8,7X.A3.11X.A3.11X,I2.12X.I2.11X.I2.15X.'2'.16X.A3.//) 
F0RMATaiX,A8.7X.A3.11X.A3.11X»I2.12X.12.11X.I2.15x. *3'.16X.A3.//) 
FORMAT(lXX.A8,7X.A3,llX.A3.11X.I2.i2X.I2.11X.I2.15X» '4',16X.A3.//) 
F0RMAT{11X.A8.7X.A3.11X.A3.11X.I2.12X.I2.11X.I2.154»'5'.16X.A3,//) 
F0RMAT(11X.A8,7X.A3,11X.A3.11X.I2.12X.I2.11X.I2.15X.'6'.16X.A3,//) 
CONTINUE 
CONTINUE 
GO TO 900 
980 V(RITE(6.695)NUMB 
697 F0RMAT(/5X. 'NUM8 = ».I4) 

GO TO 900 

55 WRITE(6.56) SEC. ROW, COL 

56 FORMAT (/5X, 'YOU HAVE MIS-SPELLED A SOIL SYMBOLl ', /5X ,' SECT ION' , 

1 I3,3X.'R0W',I3.4X,'C0L* .13) 

GO TO 51 
900 STOP 
END 


PEP03130 

PEP03140 

PEP031SO 

PEP03160 

PEP03170 

PEP03180 

PEP03190 

PEPO32O0 

PEP03210 

PEP03220 

PEP03230 

PEP03240 

PEP03250 

PEP03260 

PEP03270 

PEP03280 

PEPO3290 

PEP03300 

PEP03310 

PEP03320 

pepo;t"“ 


III 


PEP03340 
PEP03350 
PEP03360 
PEP03370 
PEP03380 
PEP03390 


ORIGINAL PAGE IS 
OF POOR QUALITlfl 


\ 



FILE. . . 

072LAKE 
AD ALIOA 
AL ALIDA 
BLABLOUNT 
BN BONO 
BP BORROW 
BR BRADY 
8S8 BREMS 
CA CARLIS 
CP CLAYPl 
DA DARROC 
DE DELREY 


LAKE 


FSL 23 
L 22 
SIL 22 
SIC 23 
XX 68 
FSL 33 
FS 34 
M 63 
C 
L 


SOIL 


68 

^ 62 
SIL 22 


DL 0ELREY5ILDCV 22 
DOA DOOR L 61 

DOB DOOR L 62 

DRB DOOR SICLS 62 
DU DUNELD XX 68 
EL ELLIOT SIL 62 
GD 6ILF0R FSL 42 
GF GILFOR MFSL 42 
GM GILFOR L 42 

GP GRAVEL XX 68 
lb LAKEBE XX 68 
LM LINWOO M 63 

LYALYDICK L 61 

LY8LYDICK L 62 

MABMARKHA SIL222 
MB MARLBO XX 66 
MH MARSH XX 68 
MM MAUMEE LFS 43 
MN MAUMEE SIL 43 
HO MILFOR SILO 22 
MR MILFOR SICL 22 
MS MILFORSICLSA 22 
MT MILFOR XX 22 
MUBHORLEY 
MUCMORLEY 
MUDMORLEY 
MUEMORLEY 
MVBHORLEY 
MVCMORLEY 
MVEMORLEY 
OAEOAKVIL 
OKBOAKVIL 
SP SANOPI 
OSAOSHTEM 
OSBOSHTEM 
OSCOSHTEM 
PC PEWAMO 
PE PEWAMOSICLCV 22 
PL8PLAINF FS 34 
PLCPLAINF 
RE RENSSL 
RN RENSSL 
RR RENSSL 
RS RENSSL 
RE RENSSL 
SPBSPARTA 


SIL 22 
SIL223 
SIL224 
SIL 26 
SICL323 
SICL324 
SICL327 
FS 38 
XX 38 
XX 68 
FSL 33 
FSL 33 
FSL 33 
SICL 22 


FS 34 
L 22 
LSS 22 
MLSS 22 
LCSV 22 
L 


_ 22 
FS 64 
SPBSPARTAFSSICL 63 
TA TAWAS M 63 

TCA TRACY L 11 

TCB TRACY L 12 

TCC TRACY L 13 

TRB TRACY LSICL 12 

TYB TYNER LFS 33 

UR URBAN XX 68 

WA WALLKI SIL 62 

WE WARNER SL 66 

WR WARNER SIL 66 

WK WATSEK LFS 64 

WL WATSEK LSMDV 62 

WO WASEON FSL 22 

WT WHITAK L 22 


PI 


OEIGWAL 
OP POOR QUALOT 



FILE. . 
PORTER 

. TWNSHP 

35N07W010213 

357 PI 

4 

ELA 

PORTER 

35N07W0102U 

4 

ELA 

PORTER 

35N07W010215 

3 

ELA 

PORTER 

35N07W010216 

4 

ELA 

PORTER 

3SN07W0I0311 

4 

ELA 

PORTER 

35N07W010411 

4 

ELA 

PORTER 

35N07W011309 

4 

ELA 

PORTER 

35N07H011310 

4 

ELA 

PORTER 

35N07W011311 

4 

ELA 

PORTER 

35N07W011312 

4 

ELA 

PORTER 

35N07W011313 

4 

ELA 

PORTER 

35N07W011316 

4 

MRS 

PORTER 

35N07W0U409 

4 

ELA 

PORTER 

35N07WOU410 


ELA 

PORTER 

35N07W01U11 

4 

ELA 

PORTER 

35N07W011412 

4 

ELA 

PORTER 

35N07W011413 

4 

ELA 

PORTER 

35N07W0114I4 

4 

PE 

PORTER 

35N07W0U415 

4 

PE 

PORTER 

35N07irf011416 

4 

PE 

PORTER 

35N07W011513 

4 

MRB2 

PORTER 

3SN07W011516 

4ALL 

PE 

PORTER 

35N07W020603 

4 

HK 

PORTER 

35N07W02070I 

3 

SO 

PORTER 

35N07W020702 

2 

SO 

PORTER 

35N07W020703 

3 

SO 

PORTER 

35N07W020704 

3 

so 

PORTER 

35N07W020801 

4 

WT 

PORTER 

35N07i*;020802 

3 

MRE2 

PORTER 

35N07W020803 

4 

HK 

PORTER 

35N07*J020901 

4 

HK 

PORTER 

35N07W020902 

3 

HK 

PORTER 

35N07W021001 

3 

HK 

PORTER 

35N07W021002 

2 

HK 

PORTER 

35N07W021010 

2 

SO 

PORTER 

35N07W021105 

1 

SO 

PORTER 

35N07W021205 

1 

SO 

PORTER 

35N07W02I206 

1 

WT 

PORTER 

35NJ07W021207 

1 

WT 

PORTER 

35N07W021209 

1 

HK 

PORTER 

35N07W021211 

1 

8AA 

PORTER 

35N07W021305 

2 

HK 

PORTER 

35N07W021308 

4 

HK 

PORTER 

35N07W021309 

4 

HK 

PORTER 

35N07W021310 

3 

BAA 

PORTER 

35N07ii/021312 

4 

BAA 

PORTER 

35N07W021315 

3 

BAA 

PORTER 

35N07W021405 

3 

HK 

PORTER 

35N07W021408 

4 

HK 

PORTER 

35N07W021410 

4 

WT 

PORTER 

35N07W021505 

1 

HK 

PORTER 

35N07W021506 

1 

WT 

PORTER 

35N07W021507 

3 

WT 

PORTER 

35N07W021605 

1 

SB 

PORTER 

35N07W021607 

iall 

SB 

PORTER 

35N07W030302 

1 

OE 

PORTER 

35N07W030303 

2 

0£ 

PORTER 

35N07W030401 

3 

MQ 

PORTER 

35N07W030402 

4 

MQ 

PORTER 

35N07W030404 

1 

WT 

PORTER 

35N07W030506 

2 

MRS2 

PORTER 

35N07W030605 

2 

MQ 

PORTER 

35N07W030606 

2 

MQ 

PORTER 

35N07W030607 

2 

WT 

PORTER 

35N07W030608 

2 

WT 

PORTER 

35N07W030T05 

2 

WT 

PORTER 

35N07W030706 

2 

WT 

PORTER 

35N07W030707 

2 

MQ 

PORTER 

35N07W030709 

3 

MQ 

PORTER 

35N07W030716 

2 

WT 

PORTER 

35N07W030806 

1 

MQ 

PORTER 

35N07rf030807 

1 

MQ 

PORTER 

35N07W030809 

3 

SB 

PORTER 

35N07rf030810 

1 

SB 

PORTER 

35N07H030811 

1 

SB 

PORTER 

35N07W030812 

1 

SB 

PORTER 

35N07W030813 

3 

SB 

PORTER 

35N07W030816 

4 

HK 



